ON POSITIVE CUBATURE RULES ON THE SIMPLEX AND 
ISOMETRIC EMBEDDINGS 



MASANORI SAWA AND YUAN XU 

Abstract. Positive cubature rules of degree 4 and 5 on the d-dimensional 
simplex are constructed and used to construct cubature rules of index 8 or de- 
gree 9 on the unit sphere. The latter ones lead to explicit isometric embedding 
among the classical Banach spaces. Among other things, our results include 
several explicit representations of (x\ + . . . + a;?)' in terms of linear forms of 
degree 2t with rational coefficients for t = 4 and 5. 



1. Introduction 

A cubature rule of degree n on the <i-dimensional simplex is a finite sum of linear 
combination of function evaluations that approximates the integral on the simplex, 
which becomes exact if the integrand is any polynomial of degree at most n. We 
are interested in constructing positive cubature rules on the simplex 

T d := {x G R d+1 : Xi > 0, . . . , x d > 0, 1 - xi - . . . - x d > 0}, 

where positive means that the coefficients in the linear combination are all positive 
numbers. Such objects are traditionally studied in the context of numerical anal- 
ysis. They are, however, closely related, even equivalent under mild restriction, to 
cubature rules on the unit sphere E> d of M d+1 , which are intimately connected to sev- 
eral other topics in mathematics, such as spherical design and isometric embedding 
of the classical finite dimensional Banach space [8| [13] . 

Our interests in this topic originally stem from a search for cubature rules on 
high dimensional sphere in the literature, which turned up only a couple of rules of 
degree 9. Since an invariant cubature rule of degree 9 on the sphere is equivalent 
to a cubature rule of degree 4 for the Chebyshev weight function on the simplex, 
we then looked for cubature rules on the simplex and found, to our surprise, no 
positive cubature rules of degree 4, albeit for the unit weight function, in higher 
dimensions [3J [TSJ Q15] . It turns out that there is a reason for this omission: most 
of the construction in the literature relies on choosing the points in advance and 
the points being chosen, devoid of other reasons, possess simple structures, such 
as centroid (^qij, • • • , ^Tf) an d points in the form (r, . . . , r, 0, . . . , 0) or in its orbit 
under the symmetric group of the simplex. However, such choices, as we shall 
prove in Section 3, do not lead to positive cubature rules of degree 4 or higher. 
Adding another point of a different nature, we then look for a practical way of 
constructing positive cubature rules. The most straightforward way of constructing 
a cubature rule of degree n is to solve the moment equations that match the integral 
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and the cubature sum for all polynomials up to degree n. A well known theorem 
of Sobolev |17j states that if the cubature rule is symmetric under certain finite 
group, then only polynomials invariant under the same group need to be considered. 
Restricting to invariant polynomials reduces the size of the equations considerably. 
The difficulty, however, lies in that the moment equations are nonlinear in the 
coordinates of the points. Such nonlinear equations can be solved mostly only by 
numerical means and for cubature rules of moderate degrees in lower dimensions; 
see, for example, [53] for the latest effort in this direction. By adding only one point 
of the form u a — (a, b, . . . , b, 0, . . . , 0), we shall provide a method of constructing 
positive cubature rules without solving the equations numerically. 

The positive cubature rules that we found are of degree 4 and 5 on the simplex 
and degree 9 and 11 on the sphere, from which cubature rules of the same degree 
but with far less number of points can be deduced from a method due to Victoir [20] 
that uses orthogonal arrays and block designs. Even after the deduction, however, 
the number of nodes for our cubature rules can still be fairly large comparing to 
non-positive cubature rules. If the positivity is not required, then there exist cuba- 
ture rules of remarkably small number of points of higher degrees for all dimensions 
O 112) . For simplex and sphere of higher dimensions, the positivity exerts a high 
demand that is not easy to meet. Besides the numerical stability that they auto- 
matically bring, the positive cubature rules are often useful in other context. The 
cubature rules of degree 2t + 1 on the sphere with central symmetry are equivalent 
to cubature rules of index 2t on the sphere, which give, when they are positive, 
isometric embeddings between classical Banach space l d l— ► ^2t- Moreover, some of 
our solutions are rational, meaning that the coefficients are all rational and a is a 
square of an integer, which lead to representations of [x\ + . . ■ + x 2 d ) t as a sum of real 
linear forms of power 2i for t = 4, 5 with rational coefficients. Such representations 
have implications in number theory. 

The paper is organized as follows. The next section is preliminary, where we 
sum up the background on cubature rules and their connection to other topics. 
The non-existence of the positive cubature rules for simple choice of points and our 
method of constructing positive cubature rules are given in Section 3. The method 
is applied to integrals with respect to the constant weight function on the simplex 
in Section 4 and to Chebyshev weight function, which connects to cubature rules 
on the sphere, in Section 5. 



2.1. Cubature rules on the simplex and on the sphere. For x <G T d , define 
\x\i = x\ + . . . + Xd- For 7 > — 1 and x G T d , define 



Let IT^ denote the space of polynomials of degree at most n in d variables. A 
cubature rule of degree n for the weight function W 7 is a finite sum of function 



2. Preliminaries 



Wj(x) = (xi . ..x d {l - |af|i)) 



x G T d . 



Let w-y be the normalization constant of W 7 defined by 
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evaluations such that 

r N 

(2.1) Wl / f{x)W 1 {x)dx = V X k f(x k ), V/ G n d , 

and there is a function /* G U^ +1 for which the above equation fails to hold. We are 
particularly interested in cubature rules that have all nodes inside T d , i.e., Xk G T d 
and that are positive, meaning that all A& > 0. 

The simplex T d is invariant under the permutations of its vertices. To capture 
the symmetry, it is often easier to work with homogeneous coordinates, £ € M. d+1 
with |£| i = 1, under which the simplex T d is equivalent to 

T d+1 := U e R d+1 : Ci > 0, . . . , £m > 0, \^ = l}. 

Under homogeneous coordinates, Wj(x) becomes £ 7 and we have 



(2.2) / f(0Cd£= / ffal-lxlJW^dx. 

JT d+1 JT d 

Working with T~ d+1 allows us look at cubature rules of a particular form. 

Definition 2.1. Let V^ 1 denote the space of homogeneous polynomials of degree 
n in R d+1 . A cubature rule is of index n (in contrast to degree n) if it is exact for 
all polynomials in V d+1 . 

Proposition 2.2. The cubature rule (|2.1[) with all nodes inside T d is equivalent to 
a cubature rule of index n for the integral against £ 7 on T d+1 : 

N 

(2.3) Wj / /(or^e = E A */^)' y f e v n +1 , 

where an d x k are related through = (xk, 1 — |#fc|i) G T d+l ■ 
Proof. Every homogeneous polynomial F G "P^ +1 can be written as 

n 

F(x 1 ,...,x d +i) = ^2fj(xi,...,x d )xj~l with fj€ll d , 
j=o 

which implies, in particular, that f(x) := F(x\, . . . , Xd, 1 — Mi) G 11^. Indeed, 
applying ()2.1[) to such / gives, by (12. 2[) . the cubature rule ()2.3|) . On the other 
hand, for each / G 11^ with < m < n, applying (|2.3|1 to the homogeneous 
polynomials f(x u x d)\x\1' m G V d+1 gives $tft. □ 

Switching from (|2.1|) to (I2.3[) has one advantage, namely, that we can disregard 
the restriction of G T d+l ■ Indeed, (|2.3[) is equivalent to 

»7 / /(or^e = E x tm/r k ), v/ g 

where A* = AfcrJ? for any > 0. Therefore we can always rescale so that the 
points are inside T d+1 ■ This fact is especially useful for constructing cubature 
rules of lower degree in high dimensions, as will be seen in the following sections. 
For a G Nq +1 , the moment with respect to £ 7 on F d+l is given by 

( j m " ' W V T ^ ^ ((7 + l)(d+l)) Wl ' 
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where (a) m =: a(a + 1) . . . (a + m — 1). 

We arc interested in two integrals on the simplex, corresponding to 7 = and 
7 = —1/2, with respect to the weight functions 



W {x) = 1 and W_ x / 2 {x) = \j\Jx x ■ -.x d (l - \x\i), 

respectively. The hrst integral is traditionally studied in numerical analysis, whereas 
the second one is of interest because a cubature rule for this weight function corre- 
sponds to a cubature rule for the surface measure a on the sphere E> d . In fact, the 
cubature rule (|2.ip of degree n is equivalent to a cubature rule of degree 2n — 1 on 
the sphere, 

„ M 

<r d f(x)da = J2fi k f(y k ), V/€ n 2 „_ 1 (S d ), 
J V fc=i 

where II„(§ d ) denotes the space of polynomials of degree at most n in d+1 variables 
and <7d = w_ 1/2/2 is the surface area of § d . To state this relation [2T] more 
precisely, for x <E R d+1 , let t(x) denote the number of non-zero components of x. 
For x = (x 1( . . . , x d ) e T d , let 



X = (sci, . . . ,Xd, 1 - \x\i) and vX = (y/xi, . . . , y/x^, \A - |x|i). 

Proposition 2.3. // £/ie cubature rule (|2.ip is 0/ degree n and it has all nodes 
x k G T , i/ien 

, d M . 

^ / /Wn^ +1 ^) = E^y E v/en 2n+1 (§ d ), 

where u k = \fX~k~. Furthermore, this relation is reversible. 

A subset X of R d+1 is said to be podal, if for any distinct x,x' £ X, x 7^ —a/. 
We denote —X = {(— xi, . . . , — x^+i) : (xi, . . . , id+i) G X}. The set X U (—X) is 
an antipodal subset of R d+1 . The following is a slight generalization of Proposition 
4.3 of [13] , which follows immediately from Propositions 12.21 and 12.31 

Corollary 2.4. Let {27, . . . , Xjv/2} ^ e a P°dal subset ofS d . If 27 and X k define a 
cubature rule of index 2t on E> d , then 

r N/2 x 

(2.5) a d / f(x)da(x) = V ^{/(x fc ) + f(-x k )}, V/ G II 2t+1 (§ d ). 

Conversely, any centrally symmetric cubature rule of degree 24+1 on S d with nodes 
±27 can be reduced to a cubature rule of index 2t on S d with nodes x k ■ 

2.2. Invariant cubature rules. Let Q be a group of bijective linear transforma- 
tions on R d+1 . Let 17 be a (/-invariant subset of R d+1 and \i be a (/-invariant measure 
jj, on fi. For a point x G R d+1 , we denote the (/-orbit of x by xQ. A function / on 
ft is Q-invariant if for any g G Q and x G SI, f(xg) = f(x). A cubature rule with 
nodes x k and weights X k , with respect to y,, is Q -invariant if x k are partitioned into 
(/-orbits, and x k g = x k i implies X k = X k > . 

We are particularly interested in the group A d of permutations of the axes of 
R d+1 . The group Ad is the symmetric group Sd+i, and an ^-invariant function 
on R d+1 is a symmetric function of 27, . . . , x d +i- We also deal with the Weyl group 
B d+ i of type B. This group contains a subgroup of reflections of the axes, which is 
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isomorphic to the elementary 2-group *Z d+1 . Throughout this paper we consistently 
use these notations. 

The weight function W 7 is invariant under any permutation of the d + 1 vertices 
of the simplex T d , fixing the origin. This means that T d is y^-invariant in the ho- 
mogeneous coordinates of T d+1 . A symmetric cubature rule for the weight function 
£ 7 := £7 • • • Cd+i i s °f the form 

r M \ 

<"> not 

where |xfc^4d | denotes the cardinality (the number of distinct elements) of the orbit 
XkAd- It is obvious that the number of nodes of (|2.6I) is Ylk=i \ x k<Ad,\- According to 
a well-known theorem of Sobolev |17j , a symmetric cubature rule of index n holds if 
and only if ()2.6j) holds for every / £ {'P d+1 ) Ad , where {V d+1 ) Ad denotes the space 
of homogeneous polynomials invariant under Ad- The dimension of this space is 
given by the partition number 

dimCP^+V* = card {/3 G N d+1 : \[3\ = d, fa > fa > ■ ■ ■ > /3 d } . 

For small n this means that we only need to verify a small number of invariant 
polynomials, which can possibly be solved by working with the moment equations. 
A basis of the space { r P d+1 ) Ad is given by symmetric polynomials of degree n, which 
are indexed by partitions. Let I = (£i, . ..,£ r ) be a partition, that is, t\ > t% > 
... > l r > 0, li £ N. We denote the symmetric homogeneous polynomial associated 
to i by 

s {e}(x) = sym^ 1 • • -x l r T ). 

2.3. Isometric embeddings of Banach spaces. Let p be a positive integer with 
p =/= 00. We define the norm | • | p : ~R d — > K by 



■r 



d N i/p 

i=l 



The Euclidean space M. d endowed with the norm | • \ p is a classical finite-dimensional 
Banach space, usually denoted by l d . 

Given two Banach spaces lp,l d ,, a classical problem in the theory of Banach 
spaces asks when there exists an R-linear map F : l d — > l d , such that 



|F(a;)|p/ = \x\ p for every x £ I 



Such a map F is called an isometric embedding from l d to l d , . To exclude trivial 
cases, we assume p' > p > 2 and p^p'. It is well known p~3j Theorem 1.1] that if 
p,p' 7^ 00 and an isometric embedding from l d to l d , exists, then p = 2 and q is an 
even integer. Hereafter we only consider the case p — 2, p' = It < 00. 

Theorem 2.5. The following are equivalent. 

(i) There exists a cubature rule of index 2t on S ld_1 with N points; 
(ii) There exists an isometric embedding l d — > l^ t ; 
(Hi) There exist N vectors ri, • • ■ , rjy £ K" 1 such that for any x £ R d , 

N 

(2.7) 5>,r 4 ) 2t Hx,z>*, 

i=l 
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where (•,•) denotes the usual inner product. 

This theorem was proved by Lyubich and Vaserstein [15] . The equivalence "(i) 
<=> (hi)" was also observed by Reznick |14] in connection with Waring's problem in 
number theory. 

Given a cubature rule of index 2t on with N nodes x^ and weights A&, we can 
explicitly construct an isometric embedding 1$ — > l^: For f(y) = (x,y) 2t g V% t , 

N 

(2.8) VyW^j-J Kx,?/)! 2 '^^) = c t (x,x)*, 

7~T Js*- 1 



7 = 1 



where 



The mapping 



c t :=a d -i y 1 da(y), y = (j/i, . . . , y d ) 



^(z) = ( 2 \/\i/c t (x,xi), 2 {/\ N /c t (x,x N )) 
defines an isometric embedding — > l% t . (|2.8|) can be transformed into the form 
(|2~7]) by setting n = 'y/Xi/ctXi. 

The value of ct can be easily evaluated (it goes back to Hilbert [5J), 

_ r(f)r(t+l) _ (2t- !)!!(<£ -2)!! 

Ct ~ r(t + f)r(i) ~ (d + 2t-2)i! ' 

which is of course a rational number. Thus, for example, if there is a /S^-invariant 
cubature rule of index 2t on S d_1 such that the weights X k are rational and the 
nodes xi, are of the form f-?^=, . . . , -4^=, -7^=, . . . , -7^=) with ak,bu, ci~ being rational, 
then (|2.8p shows that 

C{(xf + . . . + a; 2 ,)' = —(a:, (ft*, . . . , b k , c k , . . . , c k )) 2t , 

k=l Q k 

which gives a representation of (x, x) as a sum of the 2t powers of linear forms with 
only rational coefficients. Such representations are of interest in number theory; 
see [U [14] for details. 

In Section 5 we will construct some examples of index eight cubature rules on 
the sphere whose nodes and weights involve only rational numbers. 

2.4. Block designs and orthogonal arrays. Victoir [3D] proposed a method 
that can be used to reduce the size of invariant cubature rules on the simplex and 
on the sphere. His approach utilizes two combinatorial objects, namely i-designs 
and orthogonal arrays. We first explain these two concepts. 

Let t, k, v, A be integers such that < t < k < v. A t-design is a system of a set 
V of v elements, called points, and a collection B of fc-element subsets of V, called 
blocks, such that every i-element subset of V appears in exactly A blocks of B. This 
is denoted by t-(v, k, A). It is well known (cf. [7 ) that if a t-(v, k, A) with b blocks 
exists, then 

An incidence matrix of a t-(v, k, A) is a v x b zero-one matrix with rows and columns 
being indexed by points and blocks, respectively, such that for x £ V, B £ B, its 
(x, B)th entry is 1 iff x £ B. 
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Let s,l,L,/i be nonnegative integers. An orthogonal array with strength s, con- 
straints I and index fj,, is an L x I matrix, such that in every s columns, each of the 
2 s ordered s-tuples of elements ±1 appears exactly fi times among L rows. This is 
often denoted by OA(L, I, 2, s); we do not put \i in the notation, because /i = L/2 S , 
by the definition of OA. 

We need a bit more notations. Let e\ = (1, 0, . . . , 0), . . . , e^+i = (0, . . . , 0, 1) be 
the standard basis vectors of R d+1 . Define 



k d+1 
i=l i=k+l 

We recall the notation r(-) as in Subsection 2.1. Victoir proved the following. 

Proposition 2.6. (i) Assume there exists a symmetric cubature rule of index t on 
T d+1 , such that 

Moreover assume that for i = 1, . . . ,£, there exists a t-design with d+1 points 
and hi blocks of size ki, with an incidence matrix Ii. Let Xi be the columns of 
aiIi + (/3i — ai)Jd+i,bi, where Jd+i,bi denotes the all-one matrix of size (d+l)xbi. 
Then, 

J 1 i=i iex, i=£+i 1 1 xexiAd 

defines a cubature rule of index t on T d+1 ■ 
(ii) Assume there exists a -invariant cubature rule of degree 2t + l on the sphere 

[ M Ai 

o-d f(x)da(x) = V —7^t V f(x). 



x£Xi 



( = 1 



Moreover assume there exist OA(\Xi\, t(x{), 2, 2t+l) with rows Xi, i = 1, . . . , M. 
Then, 

M 



a d f f( x )da(x)=<TwiY, 



i=i 1 1 xeXi 
defines a cubature rule of degree 24 + 1 on E> d . 

Let us briefly explain on how Proposition 12.61 is used; see [501 Subsection 4.4] . 
We start with a symmetric cubature rule of index t on T d+1 ■ By Proposition l2.6l fi). 
if there is a t-design for the set {vj?"^ : i = 1, ...,£}, we can reduce the size of 
this cubature rule, where (I2.9[) is used to calculate the values of bi with u = d+ 1. 
By Propositions 12.31 and 12.21 the resulted cubature rule can be transformed to a 
Z2 +1 -invariant cubature rule of degree 2t + 1 on § d . Finally, by Proposition l2.6l fii'). 
we reduce the size of the last cubature rule without reducing its degree of exactness. 

We emphasize that the process mentioned above is not exactly the same as 
that of Victoir, who wrote his method in terms of cubature rules of degree-type. 
We slightly modify his approach to fit the cubature rules of index type, where 
Proposition 12.21 serves as a bridge between cubature of degree-type and cubature 
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of index-type. In Sections 4 and 5, we will apply the method to derive several new 
cubature rules on the simplex and on the sphere. 

3. Symmetric cubature rules for the simplex 

In the first subsection, we prove a result that explains why many of the cubature 
rules in [18] are not positive. We then present our method of constructing cubature 
rules in the following three subsections. 

3.1. Positive cubature rules. As mentioned in the introduction, many known cu- 
bature rules on T d are invariant and constructed by choosing their nodes among the 
centroid (-7==, . . . , ^j=^ ) and orbits of the points of the form (r, . . . , r, 0, . . . , 0) <G 

T d . Let us define these points more precisely. For 1 < i < d, define v k = 
e\ + ■ ■ ■ + &k — (1, • • • , 1, 0, . . . , 0). Then the points are 

(3.1) (-^=,...,-^=) and (r k v k )A d , l<k<d+l, 

where r k are real numbers such that r k Vk = (ffe, . . . , r k , 0, . . . , 0) £ T d , and {rkVk)Ad 
means the orbit of (rkVk, 1 — |?"fcUfc|) under the symmetric group of d+ 1 elements. 

Theorem 3.1. If a symmetric cubature rule of degree n with n > 4 for W 7 on T d 
uses only points form those in (|3.1[) . then it cannot be positive. 



Proof. The key ingredient of the proof is an orthogonal polynomial of degree 4. 
By Proposition 12.21 we can consider a homogeneous orthogonal polynomial, call 
it P4, of degree 4 on T d+1 ■ To derive this polynomial, we use a relation between 
orthogonal polynomials on the simplex and those on the sphere |21j . which shows 
that Y (x) :— P±{x\ , . . . , x 2 d+1 ) is an orthogonal polynomial on § d with respect to the 
measure nf=i |^i| 27+1 ^ CT - The polynomial Y then satisfies the equation A^Y = 0, 
where A/, is the Dunkl Laplacian associated with the group 1 d+1 . The operator A h 
is usually a differential-difference operator, but it becomes a differential operator 
when it is applied to polynomials that are invariant under Zj" 1 " 1 (see, [SJ p. 156]). 
Consequently, the polynomial Y is a solution of 

d+1 a2y d + 1 1 J)v 

i=i 1 fc=i 
With this explicit expression of A^, a straightforward computation shows that 

A h S {8} (x) = 16(j + 4)S {6} (x), 
&hS {6 , 2} (x) = 12(7 + 3)S {4l 3}(x) + 4d( 7 + l)S {6} (a;), 
A h S {m (x)=8{'y + 2)S {4i2} {x), 

from which it is easy to verify that 

A . o x/,Q In 4 (7 + 4) 6( 7 + 3)( 7 + 4) 

A h [S {8} + aS {6 , 2} +bS {m ]=0, a = -^_ Ty , 6= d(7 + 1)(7 + 2) . 

Consequently, as discussed above, the polynomial 

Pa = S{4} + aS 1 ^!} + bS{ 2 ,2} 
is an orthogonal polynomial of degree 4 with respect to £ 7 on T d ■ 
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A quick computation using Lemma 13.21 below shows that, for 1 < k < d, 

x£v k A d 

fc(d( 7 + 1)(7 + 2) + (7 + 4)( 7 + 17)(fc ^1)) 
d( 7 + l)(7 + 2) 

Consequently, should an invariant cubature rule of index 4 that uses only the orbits 
of Vk be positive, we would have 

r M a 

which is an obvious contradiction. Since a cubature rule of index n, n > 4, is 
automatically a cubature rule of degree 4 as seen by Proposition l2.21 this completes 
the proof. □ 



Theorem 13.11 explains why there are no positive cubature rules of degree 4 or 
higher in higher dimensions among those given in [18 . When 7 = —1/2, this result 
was proved by Bajnok [2J, where he worked with the cubature rules on the sphere. 

Accordingly, in order to construct positive cubature rules, we need to include at 
least one point not among those in (|3.1j) . We shall add one point in the form of 

ae\ + be 2 + . • • + be p+1 = (a, b, . . . , b, 0, . . . ,0), a, b > 0, 

which has p elements of b and one a. By Proposition ^. 21 wc work with homogeneous 
polynomials for the integral over T d+1 , for which we can rescale of the point and 
consider (a, 1, . . . , 1, 0, . . . , 0) with a > and a^l. 

Let ki, . . . , k q be four distinct integers, 1 < fci, . . . , k q < d + 1. We consider 
cubature rules of index t in the form 

^L/w=e4 e fM+Tn e /(»). 

3=1 V kj I x£v k .Ad x£u aiP Ad 

where u a . P := aei + ei + ■ ■ ■ + £p+i = (a, 1, . . . , 1, 0, . . . , 0) with a > and o ^ 1. 
Since the cubature rule is of index type, it can be rewritten, upon rescale, as 



(3.3) 



Vt-1 \ " r 

d+1 ^ . 



\ a + P 

which has all nodes on f d+1 and has the weights given by 

(3.4) Xj = Ajkj, j = l,...,q, and A g+ i = A q+ \ (a + pf , 

By Proposition l2.21 the cubature rule (|3.3j) is equivalent to a cubature rule of degree 
t on T d . 
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3.2. Cubature rule of degree 4. We first describe our method of constructing 
positive cubature rules for rules of degree 4 on T d . In this case we choose q = 4, 
the reason will be clear by the end of this subsection. 

Since the cubature rule (|3.2[) is invariant under the group Ad, we only need to 
work with symmetric polynomials. There are five symmetric homogeneous polyno- 
mials of degree 4 in M, d+1 for d > 3, which are 

S{4}i ^{3,1}' £{2,2}, S{2,1,1}> ^{1,1,1,1}- 

We will need the values of these polynomials on the nodes of the cubature rule p.2[) 
which come down to a straightforward counting and the results are collected in the 
lemma below. 

Lemma 3.2. Let £ = (£i,...,£ r ) be a partition oft such that iy+i = ... = 
£i j+1 ,tij 7^ where j = 0, s—1. For l<k< d+l andu cltP = (a, 1, . . . , 1, 0, . . . 

we have 



S{i}{v k ) 



r 

\ r J \}i ~ *o,«2 — h,-" - i s -i 

s-l 



r- 1 \f p 

h-io,--- -ij - i s -\) \r - I 



3=0 

+ f r wp 

\*i - *o>*2 - *ij • * ' . l s - i s -ij \r 

where (j 1 _j 0i j 2 _j 1) ... i ,- s _j 3 _ 1 ) — (ii— io)!(»2— ujl---(is— ' 

We now apply the cubature rule (|3.2[) to the symmetric polynomials. Let m^y 
denote the integral of Srn on T d+1 , 

m {e} := w y / S {e} {£)Cd£. 

JT d+1 

Their values can be easily evaluated using the explicit formula of the moments of 
£ 7 <i£; see the following section. 

Proposition 3.3. A positive cubature rule (|3.2p of index 4 exists if 

(3.5) A 5 = (m {3jl} - 2m {2i2} ) / (a(a - l) 2 p), 

and, for distinct positive integers ki,k2,k3,ki, the following system of equations 

4 

Ajkj + A 5 (a i + p) — m{ 4 } 

i=i 



£ ^ ft) + A * ( a2p + (?) ) = m ^ 2 } 



(3.6) 



i=i 

4 



3 E A i(3j +A » (V + ^Q +3(3)) =»W,i} 
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can be solved with positive Aj and a. Furthermore, if {A\, A2, A3, A4, a} is such a 
solution, then the positive cubature rule of index 4 with all nodes in J~ d+1 takes the 
form (|3.3[) with weights given by (|3.4p . 

Proof. By Sobolev' theorem, the existence of the cubature rule (|3.2[) is equivalent 
to the solvability of the system of equations 

4 

22 A 3^{l}( v k j ) + A 5 S{e}(u a , p ) = m {e} 
j=i 

for all partitions of 4. The equations in (|3.6p correspond to I = {4}, I — {2,2}, 
£ = {2, 1, 1} and £ = {1, 1, 1, 1} by Lemma [3^1 The equation for £ = {3, 1} takes 
the form 

2) ' (i a3 + a )P + 2 (^)) =^{3,1}- 

Taking the difference of this equation with 2 times of the equation for A = (2,2), 
we obtain 

a(a - l) 2 pA b = m {3A} - 2m {2 , 2 }, 
which gives the formula for A§, □ 

We should comment on the nature of this proposition. In order to find cubature 
rules that use small number of points, we should solve the moment equations with 
a, in u a , p , as a variable, which means that we should solve, heuristically, the mo- 
ment equations with five parameters, say A\, A2, A3, A4, a. However, the equations 
in these parameters are nonlinear and can only be solved numerically. Instead, 
we have added one more coefficient and solve the linear system of equations for 
Ai, A2, A3, A4, A$ symbolically with a as a free parameter. We then look for pos- 
itive cubature rules by choosing a judiciously. Whenever possible, we look for 
solutions for which Ai are all nonnegative and with as many being zero as possible 
to reduce the number of nodes. In fact, in many of our examples in the follow- 
ing section, we have at least one zero coefficient. For those, in essence, we have 
solutions for the non- linear system of five equations with five parameters. 



3.3. Cubature of degree 5. With the same set of nodes, we now consider cuba- 
ture rules of degree 5. It turns out that they are easier to deal with. There are 
seven symmetric homogeneous polynomials of degree 5 in which are 

£{5}> "5{4,1}> £{3,1,1}' "{3.2}) ^{2,2,1}! £{2,1,1,1}! £{1,1,1,1,1}- 

Their values on the potential nodes are summarized in Lemma 3.2. 

It turns out that we can solve a explicitly and turn the nonlinear system equations 
of moments into linear system. In other words, the seven moment equations can 
be solved with 7 parameters, a, A\, . . . , A§. Hence, we choose q = 6 in (|3.2[) . 

Proposition 3.4. A positive cubature rule (|3.2j) of index 5 exists if q = 6, 

f17\ „ P Z 1 Z m < 3 - 2 > 1 A m 14,l} ~ m {3,2} 

(6.i) a= — i, a 6 — — — — — — — , 

2 W{3,i,i} - m {2,2,i} pa(a-iy{a+l) 
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and, for distinct positive integers k\, fe, £3, k±, ks, the following system of equations 

5 

™{5} 



J2^k 3 +A 6 (a 5 +p) 



3=1 

5 



2 A i [o) + Ae ( a4 P + a P + 2 ( o ) ) = to ki> 



3=1 

5 



(3-8) 3 £ A, ; Q + A 6 ( (a 3 + 2a) ( J ) + 3 ( J ) ) = m {3 , 1)1} 



3= 
5 



4 X>00 + ((«' + 3 «)(s) +"(4) ) = 



3=1 



can 6e solved with positive Aj. Furthermore, if {A\, A2, A3, A4, A§} is such a so- 
lution, then the positive cubature rule of index 5 with all nodes in f d+1 takes the 
form (13.31) with weights given by (13.4[) . 

Proof. Again, by Sobolev' theorem, we need only to solve the moment equaitons 

4 

3=1 

for all partitions of 5. The equations in (13.81) correspond to I = {5}, I = {4, 1}, 
£ = {3, 1, 1}, £ = {2, 1, 1, 1} and I = {1, 1, 1, 1, 1} by LemmaO The equations for 
£ = {3, 2} and t = {2, 2, 1} take the form 

2 il A ^i) + ^ (V +« 2 )P + 2Q) = ™{3,2}, 



Taking the differences of the equations for the pair m^ij, m^^} and the equations 
for the pair of mxa i n and m.{2,2,i}j we see that 

- I? A, 



pa(a - l) 2 (a + 1)A 6 = m^,i} - ™{3,2}, (jf) a ( a ~ X ) 



i6 = "1(3,1,1} - TO {2,2,1}, 

from which we can solve for a and Ar. □ 



According to the proposition, the value of a is determined before hand so that the 
remaining moment equations become a linear system, which can be easily solved. 
The fact that we no longer have the freedom of choosing the value of a reduces our 
chance of finding positive cubature rules, but it does make our job of constructing 
cubature rules easier. 
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3.4. Cubature rule of degree 6 and above. In this subsection, we show that 
the cases t = 4, 5 are the most that we can do with the cubature rule (13.211. 



Proposition 3.5. If a cubature rule (13. 2[) of index 6 exists, then 
,„ Q s _ 2(3m {2 2;2} - W{3,2,i> + ™{4,i,i}) 

"^{3,2,1} - 0771(2,2,2} 

and it has to satisfy the quadratic equation 
(3.10) 

(™{4,2} - m {33} )a 2 + (2m {4 2} - 2m { 3 !3} - m {5A} )a + rn/4,2} - 2tt7 {3 3} = 0. 
Proof. We substitute S;s t n(x), S^ 2 }{x), <S{3,3}(x) into (|3.2I) . Then by Lemma I3~2l 

2 E 4» ( 2 ) + ( (fl5 + a)p + 2 (2) ) = TO{5 ' 1} 



9 



l{4,2} 



'/ 



E ^ ( 2 ) + ( fl3p + (l ) ) = ™{ 3 ' 3 > 



By multiplying both sides of the third equation by 2 and subtracting this new 
equation from the first equation, we obtain 

(3.11) pA q+1 a(a - l) 2 (a + l) 2 = m {5:1} - 2m {3y3 y 
Doing the same thing for the second and third equations, we have 

(3.12) pA q+1 a(a - l) 2 a = rn{4,2} - 2to{ 3i3 }. 

(|3~T0l) follows by (j3~TT]) and (|37l2jt . Next we shall show flU]). We substitute 
S{4 : i,i}(x), S{ 3 , 2 ,i}(x), S{2,2,2}( x ) into Pl2]l . Then by Lemma E21 

3 I>(Y) +4*1 (V + 2«)(f) +3^)) =m {4 , 1)1} 



3 J VT± V v — y V2 7 13 

k A , „ A,/_3 , _2 , nM , e /V 



3=1 

<I 



6 £ ( I ) + Ah-i ( 5(» 3 + " 3 + «) ( 2 J + 6 ( 3 ) ) = m {w} 



j=l v / \ \ / \ 

By multiplying both sides of the first equation by 2, and then subtracting this 
equation from the second equation, we have 

(3.13) - 2A q+1 a(a - l) 2 (a + 1) = m{ 3 ,2,i} - 2m {4il| i}. 

As for the second and third equations, a similar consideration shows 
(3-14) 2A q+1 a(a - l) 2 ( ] = m{ 3j2 ,i} - 6to {2 , 2 ,2}- 



x 2 / 

(|3~TU1) thus follows by (j3~T3l and (l3~T4]) . □ 
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In the two cases that we are interested in, 7 = and 7 = —1/2, Proposition 13.51 
can be used to show the nonexistence of cubature rules of index 6 on T d+l ■ Indeed, 
in the case 7 = 0, (|3.10j) becomes, as can be verified using (12.41) . a 2 + a + 1 = 0, for 
which there is no real solution. In the case of 7 = —1/2, (|3.9I) shows that a = 7/3 
and (13. 10[) becomes a 2 — 6a + 1 = 0, which does not have a = 7/3 as a solution. 
We sum this up in the following theorem. 

Theorem 3.6. For 7 = and —1/2, no matter what q and a are, there does not 
exist a cubature rule of the form (|3.2j) of index 6 or more. 



4. Positive cubature rules of degree 4 and 5 on the simplex 
In this section, we consider cubature rules of degree 4 for the integral 

(4.1) 1 / = i / f( x )dx. 

d\ J T d+i d\ J T d 

The moments for this integral is given by (|2.4[) with 7 = (0, . . . , 0), 

(4-2) ^StTT^ 

4.1. Cubature rules of degree 4. The integrals of for partitions of 4 arc 
given by 

24 , . 4 /d+1 



m ^ = (dTT)I (d + 1) ' TO{2,2} = U+rji V 2 

2 fd+l\ 1 /d + 1 



m {2,i,i} = 7T-TT- , . n* { i,i,i,i } = 



(d + i) 4 v 3 y (rf+i) 4 v 4 

In order to have smaller number of points for a fixed dimension, we shall choose 
p = d, that is, u a<p = (a, 1, . . . , 1). 

We report our results in cases depending on the values of (ki, fo, A3, k^). In order 
to have cubature rules with smaller numbers of nodes, we shall choose k\ = d + 1 
whenever possible and choose k\, k%, £3, k^ as small as possible otherwise. Most of 
our positive cubature rules are of the type (d + 1,1,2, m) (except in Subsection 
4.1.3), which are of the form 



IT' 

A3 ^ 



e(i,o,...,o)Ai 



E /w + t^v E /(*) 



; (ii,o,...,o)^ ^™ ^ e (± ) ...,i i0 ,...,o)^ 



2 ; /1 1 

16 



+ 5TT /W. 

:ce (a+d'a+d : ---'a+d)- A£i 

In presenting our examples, we will specify the value of a but do not give the values 
of Xi in general, since these values can be obtained by solving the linear system of 
equations (|3.6|) once a is known. We discuss the first case in some details to show 
how we arrive at our positive cubature rules. 
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4.1.1. Case (fci, fe, &3, k±) = (d + 1,1,2,3). In this case, the system (|3.6|) has the 
solution 

_ 6 -7a -2a 2 + a 3 -2d _ d(2 - (a - 2)(d - 4)) 

1 ~ a(a-l) 2 {d+l)B(d) ' 3 ~ 2a~B(d) ' 

-4a 2 + a(38-7d + d 2 )-2(6-5d + d 2 ) _ (a - d)d(d - 1) 

2_ 2a5(d) ' 4_ 6aB(d) 

where B(d) = (d + 2){d + 3)(d + 4), with A 5 given in ([53]) . For A 4 to be positive, 
we need a > 2. The positivity of A 3 holds if d = 3, 4 and for a < 2 + 2/{d — 4) if 
d > 4. The nominator of ^! is decreasing in d and it is not difficult to check that 
it is negative if d = 6 and 2<a<2 + 2/(6 — 4) =3. Thus, we can have positive 
solutions only in the case of d = 3, 4, 5. We give in each case the values of a and 
A j , where A j and Aj are related by, see (|3.4p , 

Ai = (d+l) 4 i4i, A 2 =A 2 , A 3 -2 4 A 3 , A 4 = 3 4 A 4 , \ 5 = (a + d) 4 A 5 , 

so that Y)j—i Aj = 1. The value of a that gives positive cubature rules is not 
unique. Whenever we can, we choose a so that one of the coefficient is zero and the 
resulting cubature rule has the smallest number of points among all choices of a. 

Case d = 3,4. The cubature rules have 19 and 31 nodes, respectively, with 
a = 1(38 -7d + d 2 + a/1252 - 372d + 93d 2 - 14d 3 + d 4 ). 



For d 7^ 3, 4, these a and Aj give cubature rules of degree 4 with non-positive 
coefficients with N = (^J 1 ) + C^ 1 ) + 2d + 3 nodes. 

Case d = 5. The cubature rule has 32 nodes. 

a = 4, A 1 = 0, A 2 = T L, A 3 = 0, A 4 = g A 5 = 

4.1.2. Case (fci, fo, fei ^4) = {d + 1, 1, 2, m), m = 4, 5, We will not give detail 

on how the solutions were found. To find the positive cubature rules, we consider 
the following two classes: one with A 2 = and the other with A 3 = 0. It turns 
out that positive cubature rules occur only in a few cases in each class. We obtain 
positive cubature rules only when 3 < d < 12 and d = 14, 15, 17. 

Class A 2 = 0. Positive cubature rules for d — 4,5, 6, 8, 9, 11, 12, 14, 17. For each d, 
the cubature rule has at most N^-d+i = ( d ^ n 1 ) + +d + 2 nodes. Its coefficients 
are given as in the above with A 2 = and 

-16 - d + d 2 + 18m - 2dm + v 7 ^ 

a = ; : With 

4(-l + m) 

A = -16(-2 + d)(d- m)(-l + m) + (16 + d - d 2 + 2(-9 + d)m) 2 . 

It is easy to see that a > 0. We obtain positive cubature rules in the following 
cases: 

(d, m) = (4, 4), (5, 4), (6, 4), (8, 4), (9, 5), (11, 6), (12, 6), (14, 7), (17, 8). 

The case d = 4 is noteworthy, since its number of nodes, 20, is less than the one 
in 4.1.1; moreover, Ai = and all coefficients are rational numbers, as seen by 
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Case d = 4, m = 4. The cubature rule has 20 points. 

2 32 25 

a = 6, Ai=0, A 2 =0, A 3 = — , A 4 = — A 5 = — . 

21 63 63 

Another case that is noteworthy is d = 12, which has both X 2 and A3 equal to 
zero. 

Case d = 12, m = 6. The cubature rule has 1731 points. 

n \ 2197 \ n i n i 99 1 729 
a = 6, Ai = — — — , A 2 = 0, A 3 = 0, A 4 = — — , A 5 = 



12250' 245' 1750 

Class A 3 = 0. Positive cubature rules for d — 7, 10, 12, 15. For each d, the cubature 
rule has at most Nj-d+i — ( _t ) + 2d + 3 nodes. Its coefficient Aj are given with 

2(d-m) 
a ~ d + 2 - 2m ' 

which is positive if m < (d + 2)/2. We obtain positive cubature rules for the 
parameters: 

(d,m) = (7,4), (10,5), (12,6), (15,6). 

4.1.3. Case (fci, k 2 , £3, fc 4 ) = (d + 1, 1, 3, m). Positive cubature rules occur for sev- 
eral pairs of (d,m). We can again consider two classes, A2 = or A3 = 0. In the 
case of A2 = 0, we obtain positive cubature rules for 

(d, m) = (8, 5), (9, 5), (11, 5), (12, 6), (13, 7), (14, 8), (16, 8). 

In the case of A3 = 0, we obtain positive cubature rules for 

(d,m) = (10,5), (12,6), (15,8). 

In addition, for the parameters (d, m) = (17,8), we found a positive cubature rule 
with Ai = 0. We shall not give these rules explicitly. 

Combining the above cases, we have found positive cubature rules of degree 4 for 
d = 3, 4, . . . , 17. On the other hand, we have tried various of other combinations of 
(fci, &2, A)3, £4), but found no positive cubature rules for d > 17. 

Although our main goal is to construct positive cubature rules, it should be noted 
that the cubature rule for the case (fci, £3, k±) = [d+ 1, 1, 2, d), non-positive as 
it is, has the number of nodes in the order of d 2 /2, which is the lowest among all 
choices of (fci, k 2 , £3, fc 4 ). Furthermore, the case can be solved so that A3 = with 
a choice of a = 2, which is the case discussed in [3]. 

4.2. Cubature rules of degree 5. We use Proposition 3.4 to hnd cubature rules 
of degree 5 for the integral (|4.1I) . By (|3.7p with 7 = 0, we see that 

6(p-l) , , 12d 

a = — 1 and Ar - 



pa(a - l) 2 (a + 1) ' 



With these parameters, the equation Q3.8P can always be solved and yield cubature 
rules of degree 5 on the simplex. There are, however, only a few rules that are 
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positive. These occur at the dimension d = 3, 4, . . . , 11. All of our positive cubaturc 
rules are of the type (d + 1, 1, 2, m), which are of the form 

1 ffm = *f + & E /(*) 



d! 



ase(l,0,...,0).Ad 



1 2 j , e (i,i,o,..,o)^ 1 3 j . e (§,|,| f o....,o)* 



— ,o....,o)^t d ^ /Vp/ a;ef-T-,-T- ,...,-r- ,o....,o)^ d 

v"fc m / \ a+p a+p a+p J 

where are normalized weight so that y^-, Ai = 1, which are given by 
Ai = {d+lfA u X 2 =A 2 , A 3 =2 5 A 3 , A 4 = 3 5 ^ 4 , A 5 =to 5 A 5 , A 6 = (a + p) 5 A 5 . 

The cubature rule can be easily translated to a cubature rule on T d by Proposition 
2.2. The parameters for our positive cubature rules of degree 5 are given by 

(d,p, m) =(3, 3, 4), (4, 4, 4), (5, 4, 4), (6, 5, 4), (7, 6, 5), (8, 6, 5), 

(9, 7, 5), (10, 8, 5), (11, 11, 6), 

and the number of nodes of the cubature rule is given by 

x^ + ii + r l+ ( d : i ) + ( d : i ) + ( d+i ) 



2 J \ 3 J \ m J \p j 

It turns out that all solutions are rational. The cases d = 3 and d = 11 are notable 
as each contains zero coefficients, which are given below: 

Case (d,p,m) = (3,3,4): 19 points with a = 5 and 

Ai = , A? — — , A 3 = — , A4 — , A p; = 0. An = . 

105' 70' 35' 350' ' 525 

Case (d,p,m) = (11,11,6): 1014 points with a = 5 and 

A -0 A - 1 A - 11 A -0 A - 891 A - 4096 
Al ~ 6825' Xs " 1365' A4 " °' As ~ 2275' Ae " 6825' 

4.3. Cubature rules with fewer nodes. In this subsection we shall apply Vic- 
toir's method to cubature rules on the simplex which are obtained in Subsections 
4.1 and 4.2. 

We give three tables to list up the resulting cubature rules; the first two for 
cubature of index 4 and the third one for cubature of index 5. In each table, 
the values of d, m in the first column are the same notations as in the previous 
two subsections. Nq-d+i in the second column denotes the number of nodes of the 
corresponding symmetric cubature rules. Nfd+i in the third column denotes the 
number of nodes of cubature rules which are obtained through Victoir's method. 

We shall choose the cubature rule with fewest nodes for each (d, m). For example, 
in the case of (d, m) = (12, 6), the cubature rule given in Class 4.1.3 has 586 nodes, 
whereas, the cubature given in Class 4.1.2 has 599 nodes. So in this case we choose 
the former cubature, and do not list up the latter cubature in the table. 

The cubature rules with Nq-d+i nodes in Table Q] arise from a 4-(12,6,4), a 
4-(13,6,12), a 4-(15,7,20), a 4-(18,8,21), a 4-(ll,5,l), a 4-(16,6,6). Similarly, 
we use a 4-(12, 5, 4), a 4-(14, 7, 20), a 4-(15, 8, 40), a 4-(17, 8, 15), a 4-(16, 8, 60) to 
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obtain cubature rules with N^-d+i nodes in Table[2] Moreover we use the 5-(12, 6, 1) 
to obtain Table [3] since nontrivial 5-designs must have at least 12 points, we apply 
Victoir's method only for the case (d,p, to) — (11, 11, 6). See [7] for the existence of 
4- and 5-designs that are used here. To determine the value of N T d+i , we calculated 
the number of blocks of each design by using (|2.9[) . 

Table 1. Reducing the size of cubature of Class 4.1.2. 



(d, to) 


N T a+i 


N T d+i 


(11,6) 


1003 


211 


(12,6) 


1730 


586 


(14,7) 


6556 


901 


(10,5) 


485 


89 


(15,6) 


8041 


761 



Table 2. Reducing the size of cubature of Class 4.1.3. 



(d, to) 




N T d+i 


(11,5) 


1025 


629 


(13,7) 


3811 


951 


(14,8) 


6906 


1251 


(16,8) 


25008 


1208 


(15,8) 


12903 


1593 


(17,8) 


43794 


954 



Table 3. Reducing the size of cubature of Class 4.2. 



(d,p,m) 


N T d+i 


N T d+i 


(11,11,6) 


1014 


222 



Example 4.1. We shall describe the 11 -dimensional cubature rule of TableUi 

(TrJJT / ^ f(€)d£ = TT§kf ((iT' • • • ' IT")) + T2M2 

r11 Ee((i,o,...,o))^. ia 



xex ^(,11 i,\ , 

a;e ((3'T5'-'T5)J- /l10 



Here X denotes the set of 66 columns of a 11x66 incidence matrix of the 4-(ll, 5, 1) 
design. Each element of X is represented as a cyclic shift of the following six types 
of vectors: 

(1,1,1,1,0,0,0,0,0,1,0), (1,1,1,0,1,0,0,1,0,0,0), (1,1,1,0,0,1,1,0,0,0,0), 
(1,1,0,1,1,0,0,0,1,0,0), (1,1,0,1,0,1,0,1,0,0,0), (1,1,0,0,1,0,1,0,0,1,0). 
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Remark 4.1. The above 4-(u, k, A) designs are chosen so that A is the minimum for 
given v, k. For example, when (v, k) = (12, 6), Eq. (|2.9|) implies A = (mod 2). The 
nonexistence of 4- (12, 6, 2) was shown by B. D. McKay and S. P. Radziszowski in 
1996. The existence of 4-(12, 6, 4) is known and hence a 4-(12, 6, 4) is the minimum 
design. Some nonexistence criteria are derived from equations of block intersection 
numbers by E. Kohler (1985) and J. Bolick (1990), and others from a linear pro- 
gramming bound by Delsarte (1973). The authors knew these informations through 
discussion with Reinhard Laue. Of course, a 5-(12,6, 1) is the minimum design. 

5. Positive symmetric cubature rules of degree 9 on the sphere 

In this section we construct cubature rules of index 8 on the sphere E> d , or degree 
9 on § d_1 , that are invariant under the group Bd+i- According to Propositions 12.31 
and 12. 21 this is equivalent to finding symmetric cubature rules on the simplex for 
the integral 

Wi [ /(g) 7= t= =Wi f f{x) '''' 



(5.i) = 2 , + i; 2 d+1 ^ 



r^ 1 vCi ' ' 2 J Td v^i ' ' ' x di l ~ Mi) 

The moments for this integral is given by (|2.4[) with 7 = (— i A), 

(—5) _ (2) "! ' ' ' (2)^ +1 

(d+l\ 
{ — )\a\ 

The integrals of S^y for partitions of 4 are given by 

105 . j , 9 fd + 1 

mw = i6(fj; (d+1) ' m ™ = wnA 2 

9 fd+l\ I fd+1 

W{2 - M} = iW)A 3 J' m{1 ^ 1} = n^{ 4 

We report those results in cases depending on the values of (fci, /C2, &3, ^4)- In 
order to have cubature rules with fewer nodes, we shall choose k\ = d+1 whenever 
possible and choose ^1,^2,^3,^4 as small as possible otherwise. The values of 
p, to, a are parameters and will be chosen so that the cubature rules are positive. 
All cubature rules are reported as cubature rules of index type on the sphere, 
because of their connection with other topics. The results that we report below are 
mostly of the type 

(fei, fc 2 , k 3 , fc 4 ) = (d + 1, 1,2, to), m — 3,4, .... 

In other word, they are of the form 

f(x)da(x) = J2 f ^ + W+-T) E I{X) 

' - - N V ' 2:6(1,0,- ,0)B d+ i 



xe ( 73+T--V3+T J 8<i+1 



+ t7^tty E m + E 

V 2 y / 1 1 \ V m / / 1 1 \ 
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where we have normalized the nodes so that they are points on the sphere and 
normalized our cubature coefficients by 

Ai = A^d+1) 4 , A 2 =A 2 , \ 3 = 2 4 A 3 , \ 4 = m 4 A i , \ 5 = (a+p) 4 A 5 , 

where Ai satisfy (|3.6[> . so that X\ =1 A& = 1> fr° m which the cubature rules of 
degree type on the sphere can be deduced easily. 

The solutions are obtained in the same manner as in Section 4. Hence, we shall 
state our results without further discussion on how they are obtained. Again, the 
value of a that gives positive cubature rules is not unique. We often choose a so that 
one of the coefficient is zero to reduce the number of points. It is also interesting 
to choose a as a square of an integer, which would lead to a rational solution that 
can be used to write (x,x) 4 as a sum of 8th powers of linear forms with rational 
coefficients, as discussed in Subsection 2.3. 

For each cubature rule, we give the value of a and the parameters d,p, m, but 
we will not give the coefficients Ai, . . . , A5, which can be determined from solving 
the linear system of equations Q3.6P using the given parameters. For cubature rules 
of the same dimension, we report only rules that have either smallest number of 
nodes or have some other feature, such as a rational solution. 

5.1. Case (fci, fc 2 , £3, £4) = (<£+ l,l,2,m) and p — d. For this case we see, after 
many trails, that the system Q3.6P has a positive solution in two classes of parame- 
ters, one with A2 = 0, the other with A3 = 0. 

5.1.1. A2 = 0. Positive cubature rules hold for 

3<d<17, and d = 19, 20, 23, 24, 28, 32, 
with the number of nodes 

Nsd = 4 ( d + ^ + 2™ ( d + ^ + 2 d +\d + 2), 

or equivalently, Nj-a+i — d + 2 + ( d ^ 1 ) + ( ) in terms of cubature rules on the 
simplex. The solution is given with 

(-43 + d 2 + 45m - 3dm + y/E) 

a = - , with 

b(m — 1) 

A = -4(-3 + 3m)(-6d + 3d 2 + 6m - 3dm) + (43 - d 2 - 45m + 3dm) 2 . 

The positive cubature holds for (d, m) in the following pairs: 

m = 3:3<d<5, m = 4:5<d<9, m = 5:9<d<14, 
m = 6:13<d<17, m = 7:d=19,20, m = 8 : d = 23, 24, 
m = 9 : d = 28, m = 10 : d = 32. 

There are no positive solutions for m > 10. 

5.1.2. A3 = 0. Positive cubature rules hold for 

d = 6, 7, 10, 14, 17, 18, 21, 22, 25, 26, 29, 33, 36, 37, 40, 44. 
Furthermore, the solutions are divided into three groups: 

d = 4m — 8, 4m — 7, 4m — 6. 
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The number of nodes of the cubature rule on the sphere is equal to 

N sd =2{d+l) + 2 m ( d+1 ) + 2 d+1 (d + 2), 
\ m J 

or equivalently, Nq-d+i — 2d + 3 + { d ^ n 1 ) in terms of cubature rules on the simplex. 

Case d = 4to — 6. We have positive cubature rules if 3 < m < 8 and a — 9. The 
parameters are 

(d, to) = (6, 3), (10, 4), (14, 5), (18, 6), (22, 7), (26, 8). 



Case d = 4m — 7. The cubature rules are positive if 6 < to < 11 and 

3(3to - 7) 

a = , 

to — 3 

The parameters are 

(d, to) = (17, 6), (21, 7), (25, 8), (29, 9), (33, 10), (37, 11). 

Case d — 4to — 8. The cubature rules are positive if to = 11, 12, 13 and 

3(3to- 8) 

a = . 

to — 4 

The parameters are 

(d,m) = (36, 11), (40, 12), (44, 13). 

5.2. Case (fci, fc2, k%, hi) = (d + 1,1, 2, to), p = d and a = 9. The cubature rules 
are positive for 3 < d < 28. The parameters are 

to = 3 : 3 < d < 6, to = 4 : 7 < d < 10, to . = 5 : 11 < d < 14, 
TO = 6:15<d<18, m = 7:19<d<22, to = 8 : 13 < d < 26. 

The cases d = 6, 10, 14, 18, 22, 26 have already appeared in Case d = Am — 6 of 
Subsection 5.1.2, in which A3 = 0. In all other cases, none of \ is zero and the 
number of nodes of the cubature rule on the sphere is equal to 

N sd =2(d+l)+4^ + 2 1 ^+ 2™ ( d + ^ + 2 d +\d + 2), 

or equivalently, N-j-d+i = 2d + 3 + + ( i 1 ) m terms of cubature rules on the 

simplex, which are more than the number of nodes for cubature rules in Subsection 
5.1. However, the fact that a = 3 2 means that we obtain rational solutions for the 
cubature rules on the sphere. Therefore, as noted in Subsection 2.3, we obtain a 
representation of (X^=i x 1 ) 4 as a sum of 8th powers of linear forms with rational 
coefficients. Those identities are not only of interesting itself, but played a key role 
in Hilbert's proof of the Waring problem. See |H p. 722] for details. 
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5.3. Case (fci, fe, &3, £4) = (d+1, 1, 2, m) and 2 < p < d— 1. The positive cubature 
rules are found mostly with p = m — 1, m, m + 1. We consider two cases. 

5.3.1. p = m — 1. Here we found positive cubature rules with 

a = — (-49 - 6d + d 2 + 51m - 3dm + \fK 
6d \ 

with A = 36d(2 + d- m)(-3 + m) + (-49 + d 2 + 51m - 3d(2 + m)) 2 , 

for which A2 = 0, and the parameters are 

(d, m) =(3, 3), (4, 3), (5, 4), (6, 4), (7, 5), (8, 5), (10, 6), (12, 7), (13, 7), 
(14, 7), (15, 8), (16, 8), (17, 9), (18, 9), (19, 9), (20, 9). 

The number of nodes of such a cubature of dimension d is 

N s ,=2^ + A( d+ o 1 ) + 2-( d+1 )+2-(d+l)( d 1 
\ 2 J \ m J \m — 1 

which however is often larger than the number of points of the corresponding cu- 
bature in Subsection 5.1. 

5.3.2. a = 4. In this case we obtain positive cubature rules for the parameters 

(d, m,p) =(3, 3, 2), (4, 3, 2), (5, 4, 3), (6, 4, 3), (7, 5, 4), (8, 4, 4), (9, 4, 5), 
(10, 4, 5), (11, 5, 6), (12, 6, 6), (13, 6, 7), (14, 6, 7), (15, 6, 7). 

The number of nodes of such a cubature of dimension d is 

N sd =2(d+l) + 2 d+1 +4( d + 1 ) +2 m ( d+1 ) +2 m (d+l)( d 1 

\ 2 J \ m J \m — 1 

The choice of a = 2 3 means that all these cubature rules have rational coefficients. 
The case of (3, 3, 2) is of particular interesting, as it corresponds to a formula of 
Hurwitz (see below) that represents (x 2 + x\ + x 2 + x 2 ) 4 as a sum of 8th powers of 
real linear forms. We give the coefficients of this case explicitly. 

Case (d,m,p) — (3,3,2). 

2 11 27 

a = 4, A 1 = -, A 2 = A 3 = A 4 = 0, A 5 = -. 

Remark 5.1. (i) Using (|2.7[) . our solution in Case (d,m,p) = (3,3,2) above leads 
to the identity 

5040(Vx 2 J =6^(2a; 4 ) 8 + 60^(x l ±a; J ) 8 

i=l 4 12 

+ ^(2Xi ± Xj ± Xk) 8 + 6 ± x 2 ± x 3 ± x^) 8 . 

48 8 

This identity was found by A. Hurwitz (cf. 4, p. 721]) in connection with the 
Waring problem. We recall a 144-node cubature rule of degree 9 on S 3 given in 
Case (m,d) = (3,3) of Subsection 5.3. This is centrally symmetric, and so can be 
reduced by half by Corollary 12.41 The resulting 72- node cubature rule of index 8 
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gives the identity 

315 ( ^(2x,) 2 J = 60 22(xi ±x 2 ±x 3 ± x 4 f + 60 J2(2xif + 6 2j(2xi ± 2 Xj f 

i=l 8 4 12 

+ ^2(2xi ± 2xj ± 2x k f + ^(3^ ± Xj ±x k ±xi) 8 1 

16 32 

which implies that the 4th power of a sum of four squares of even integers, multiplied 
by 315, is a sum of at most 840 8th powers of integers. Hurwitz's identity implies 
that the 4th power of a sum of four squares, multiplied by 5040, is a sum of at most 
840 8th powers of integers. Without requiring rationality, there are some cubature 
rules with fewer nodes. Salihov [13] proved that the 120 vertices of the regular 
poly tope in R 4 give a cubature rule of degree 11 on § 3 . By Corollary 12 .4[ this is 
reduced to a cubature of index 10 with 60 nodes, which can be viewed as an index- 
eight cubature rule. On the other hand, the 4-dimensional cubature rule of degree 
9 given in Subsection 5.1 has 136 nodes. This cubature is centrally symmetric and 
so can be reduced to a 67-point cubature rule of index 8. We do not know other 
examples of index-eight cubature rules on S 3 with fewer nodes than Hurwitz's rule. 
Few explicit constructions of spherical cubature rules with few nodes are known 
even for small d. 

(ii) Shatalov [16l p. 179] found a cubature rule of index 8 on § 4 with 300 nodes. 
To obtain this cubature rule, she developed a recursive construction [161 Theo- 
rem 4.4.9] to push up the dimension of a given cubature rule. The 300-node 
cubature rule comes from Salihov's cubature rule on § 3 with 60 nodes. The 5- 
dimensional cubature rule given in Subsection 5.1 has 312 nodes, which is however 
centrally symmetric and, therefore, can be reduced to a 156-node cubature rule by 
Corollary 12.41 Thus, our cubature rule has fewer nodes than Shatalov's cubature 
rule. By Corollary 12.41 the cubature rule on S 4 given in Subsection 5.2 can be 
reduced a 161-node cubature rule, which is also smaller than Shatalov's cubature. 
This cubature rule also has rational coefficients and nodes and gives the identity 

161280(J> 2 ) 

i=l 

= 59 ^2(xi ± x 2 ± x 3 ± Xi ± x b f + 52 J2(2xi) 8 + 16 ^(2^ ± 2 Xj f 

16 5 20 

+ 4^(2xj ± 2xj ± 2x k f + 28561 ^(3xi ± x 2 ± x 3 ± x 4 ± x 5 ) 8 . 

40 80 

We also note that the 5-dimensional cubature rules given in Subsection 5.3 can 
be reduced to a 196-node and a 201-node cubature rule, which are smaller than 
Shatalov's rule. 

(iii) Lyubich and Vaserstein [13j Theorem 2.8] proved that there exists a cubature 
rule of index 2t on § rf with at most dimP 2 d t +1 (S rf ) = ( rf + 2t ) nodes. Applying Victoir's 
criterion to the cubature rules given in this section produces cubature rules with 
fewer nodes than the cubature of Lyubich and Vaserstein. For example, let us 
consider the case d — 14. It is well known (cf. [11] Table 12.1]) that there is an 
OA(2 12 , 14, 2, 8), say A. This OA has antipodality, namely, if x is a row of A, then 
so is —x. This is because the OA is constructed from a cyclic linear code over F2 by 
replacing 0, 1 e F 2 by ±1. We can easily construct an 0^4(2 13 , 15, 2, 8) by defining 
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the 2 13 x 15 array to be 

A 1 

A -1 

This OA is again antipodal. The existence of a simple 4-(15, 5, 2) design, namely, 
a 4-(15,5,2) design without repeated blocks, is known. This information is not 
covered in [7]. The authors knew this information through private communication 
with Reinhard Laue. By ()2.9[) . the number of blocks is equal to 546. Thus, applying 
Proposition [U to the OA{2 13 , 15, 2, 8) and 4-(15,5,2) design, the 620414-node 
cubature given in Class 5.1.2 can be reduced to a cubature rule of degree 9 on 
S 14 with 148574 nodes, or a cubature rule of index 8 on § 14 with 74287 nodes. 
When d = 14, the cubature of Lyubich and Vaserstein has 319770 nodes. Other 
15-dimcnsional spherical cubature rules can be constructed by applying Shatalov's 
recursive construction to a 14-dimensional cubature rule with 89964 nodes found 
by de la Harpe et al. [10]. The nodes of Harpe's cubature rule arise from shells of 
the Quebbemann lattice, and therefore are antipodal. Hence by Corollarv l2.4[ this 
cubature rule can be reduced to a half-node cubature, to which we apply Shatalov's 
recursion. The resulting cubature has 224910 nodes, which is bigger than our 74287- 
node cubature rule. We do not known other explicit constructions of cubature rules 
of index 8 on S 14 with at most 74287 nodes. 

It would be interesting to study in detail the geometric configuration of points 
of cubature given in this section, and compare with cubature obtained from the 
lattice approach. This will be discussed in a forthcoming paper. 

6. Positive symmetric cubature rules of degree 11 on the sphere 
Using the moments given in (|5.ip . the quantities of a and Aq in (|3.7p are given 

by 

(6.1) a = 1054-1, A & = 60d 



d-1 " pa(a- l) 2 (a + l)(d + 9)(d + 7)(d + 5)(d + 3)' 

The results that we report below are mostly of the type 

(h,k 2 ,k3,ki,k 5 ) = (d + 1, 1, 2, 3, m), m — 3,4, .... 
In other word, they are of the form 

**js = E /(*> + 2(£i) £ /(«) 

& l l \„ v ' ze(i,o,-,o)B d+1 
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where we have normalized the nodes so that they are points on the sphere and 
normalized our cubature coefficients by 

Ai = A 1 {d + l) 5 , A 2 = A 2 , A 3 = 2 5 A 3 , A 4 = 3 5 v4 4 , A 5 = to 5 A 5 , A 5 = (a + p) 5 ,4 6 

and Aj satisfy (|3.8p . so that J2k=i Afc = 1, from which the cubature rules of degree 
type on the sphere can be deduced easily. 

The case of d = 2p — 1 is of particular interest, since in this case a = 4 and the 
solution is rational. 

Case (d+l, 1, 2, 3, to). We found positive cubature rules for 3 < d < 23. We specify 
the parameters d,p,m. The solutions turned out to be all around d — 2p— 1. Those 
with odd dimensions are 

(d,p, to) =(3, 2, 4), (5, 3, 4), (7, 4, 4), (9, 5, 5), (11, 6, 5), (13, 7, 5), 

(15. 8. 5) , (17, 9, 5), (19, 10, 5), (21, 11, 5), (23, 12, 4), 
which are all rational solutions. Those with even dimensions are 

(d,p, to) =(4, 3, 4), (6, 4, 4), (8, 5, 5), (10, 6, 5), (12, 7, 5), 

(14. 8. 6) , (16, 9, 6), (18, 10, 6), (20, 11, 6), (22, 12, 6). 

Among these solutions, we single out the following ones because they have at 
least one coefficient zero. In each case we give the coefficients explicitly. 

Case (d,p,m) = (3,2,4). 

2 11 27 

Ai = — , A 2 = — , A 3 — -, A4 = 0, A 5 — 0, A 6 — — ■ 
15 15 8 40 

Case (d,p,m) = (7,4,4). 

A - 32 A - 1 A - 2 A - A - 2 A - 32 
Al "315' A2 -105' A3 "45' A4 ~°' A5 ~15' A6 ~45' 

Case (d,p,m) = (21, 11,5). 

14641 61 405 2625 _ 16875 

Ai — , A9 — . A 1 } — U, A4. - — , Ac; — . Af; — . 

705432' 192192' A ' 31616 1414415 18304 

Case (d,p,m) = (23,12,4). This is a particularly interesting case, since three 
coefficients are zero. 

1 506 11776 

Ai — (J, A9 — , As — (J, A4 — 0, As — , Ar — . 

4095' 12285' 12285 

Remark 6.1. (i) By (|2.7[) . our solution in the case (d,m,p) = (3,2,4) above leads 
to the identity 

22680(^a;f) 5 = 9^(2a; 4 ) 10 + 180 ^(a; 4 ± a^) 10 + J^(2x t ± ± x k ) w 

4 4 12 48 

+ J2(2 Xl ± a;j ± x fc ) 10 + 9 ^(xi ± x 2 ± x 3 ± x 4 ) 10 . 

48 8 

This identity was found by I. Schur [4, p. 721] in connection with the Waring 
problem. It is interesting to note that the coefficients of the cubature rule of the 
case (d,p,m) = (3,2,4) are exactly the same as the case (d,p,m) — (3,3,2) in 
Section 6. Namely, Schur's identity implies Hurwitz's identity, though they are 
different from each other as representations; see Remark 15.11 (i). Schur's identity 
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was given in a paper by Landau in 1908 who did not mention its connection to 
Hurwitz's identity. 

(ii) Harpe et al. [TU] constructed several even-dimensional cubature rules of index 
10 on spheres, using extremal ^-modular lattices. Under some conditions, shells 
of such lattices naturally give spherical cubature rules with special degrees. For 
example, according to Bachoc-Venkov [TJ Corollat 4.1], shells of extremal 3-modular 
lattices in M. d form a cubature rule of degree at least 5 if d = 0, 2 (mod 12). Harpe 
et al. did not discuss the explicit constructions of cubature rules of index 10 in 
odd dimensions, as well as in some even dimensions up to 24. Whereas, in this 
section, we directly constructed several cubature rules in those dimensional spaces 
(c.x. M 18 ). Of course, the lattice approach may still work well in those dimensional 
spaces. For example, collecting suitable shells of an extremal 3-modular lattice in 
R (such a lattice actually exists), we may still be able to obtain cubature rules of 
index 10 on § 17 . 
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